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We consider the gravitational lensing produced by long cosmic strings formed in a GUT scale 
phase transition. We derive a formula for the deflection of photons which pass near the strings 
that reduces to an integral over the light cone projection of the string configuration plus constant 
terms which are not important for lensing. Our strings are produced by performing numerical 
simulations of cosmic string networks in flat, Minkowski space ignoring the effects of cosmological 
expansion. These strings have more small scale structure than those from an expanding universe 
simulation — fractal dimension 1.3 for Minkowski versus 1.1 for expanding — but share the same 
■ qualitative features. Lensing simulations show that for both point-like and extended objects, strings 

| produce patterns unlike more traditional lenses, and, in particluar, the kinks in strings tend to 

. generate demagnified images which reside close to the string. Thus lensing acts as a probe of the 

small scale structure of a string. Estimates of lensing probablity suggest that for string energy 
densities consistant with string seeded structure formation, on the order of tens of string lenses 
should be observed in the Sloan Digital Sky Survey quasar catalog. We propose a search strategy 
in which string lenses would be identified in the SDSS quasar survey, and the string nature of the 
lens can be confirmed by the observation of nearby high redshift galaxies which are also be lensed 
00 ' by the string. 
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The current understanding of the formation of large scale structure presumes that clusters, walls, filaments, etc. 
\ result from the gravitational collapse of small deviations from homogeneity in the early universe, but the source 
0^ ' of these fluctuations is, at present, still a mystery. The most popular model is inflation, where vacuum energy in 
the early universe drives a period of exponential growth which stretches small quantum mechanical fluctuations to 
macroscopic scales much greater than the horizon at that time. As the universe expands, these fluctuations fall back 
into the causal horizon and collapse into the structure we see today. However, there is yet no definitive evidence in 
either accelerator experiments or cosmological observations to prove that a field with the appropriate characteristics 
exists to induce such an inflationary epoch. It is worthwhile, then, to consider alternatives to inflation which can also 
explain the origin of structure, specifically, topological defects. During cosmological phase transitions, topological 
L* ■ defects can form where adjacent fields take on different vacuum states which can only be continuously connected if a 
region of false vacuum is trapped between them. The energy trapped in these defects can then gravitationaly induce 
perturbations which will collapse into the observed structure. A particularly interesting subclass of defect models is 
cosmic strings, which are lineal structures formed when a U(l) symmetry is spontaneously broken. If these strings 
were formed at the GUT scale, they can have sufficient energy, and therfore mass, to drive the necessary perturbations 
in the ordinary matter which collapse into the structure we observe today. The question then becomes, how can we 
observe strings and compare these models with inflation? Inflationary models may best be tested by microwave 
background observations, as these models make precise predictions for the spectrum of microwave perturbations from 
180° scales down to several orders of magnitude smaller. For string models, things are less certain. On large scales, 
strings predict a spectrum similar enough to inflation as to be statistically indistinguishable jij. On smaller scales, 
issues like decoherence which smooth the acoustic peaks, are not fully understood, and while it is likely that we 
could rule out inflation if strings are the true model, we may not be able to confirm the existence of strings through 
microwave observations, given the current uncertianty about the CMBR fluctuations they produce. In this paper, we 
shall discuss another method by which cosmic strings can be observed, namely through gravitational lensing. Unlike 
other defect models, where we expect perhaps only one defect to remain in a horizon volume, simulations suggest 
that a significant length of string is observable today. The gravitational fluctuations these strings induce will bend 
light, making it possible to observe a cosmic string if it is back lit by a visible source. Since there are many sources, 
including galaxies and quasars, which can light up a string, we will consider the structure of images that would arise 
in such lensing systems containing long cosmic string. Our estimates of the string lensing probability suggest that 
string lenses occur at a rate of 10% to 30% of that of galaxy lenses in the case of quasar sources, a few dozen of which 
have been observed. We show that the resulting images have a unique signature, that is, along with the image pair 
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that is associated with infinite straight strings, there will usually be a series of smaller, demagnified images which 
reside closer to the string itself. We conclude that new quasar surveys with large angular sky coverage, like the Sloan 
Digital Sky Survey, should contain a significant number of string lenses and be able to definitively observe or rule out 
the cosmic string model, at least for flat space. 

This paper is organized into several sections. In the following, we estimate the lensing probability for long cosmic 
strings and compare it to that for galaxies. In §111 we calculate the deflection of a photon in the presence of a long 
cosmic strin g. In § |Tv| , we discuss the basic theory of gravitational lenses. The next sec tion is broken into two parts, 
where § V A discusses the generation of the long strings used in our calculations and § V B explains how these were 
used to find lensed images. Finally, in §VI we discuss our results and suggest a search strategy for finding string 
lenses. 



II. LONG STRING LENSING PROBABILITY 



We begin by considering the likelihood observing a string lens system. Quasars represent the best objects for seeing 
long string lensing, as they are bright, with images not likely to be lost in any background, and new surveys like 
the Sloan Digital Sky Survey will observe them in large numbers (on the order of 10 5 quasars over 1/4 of the sky). 
We would like to know how many quasar-string lensing systems we can expect to observe in such a survey, but this 
depends on how a particular quasar sample is observed. A simpler calculation is to measure the optical depth for 
lensing, that is the the probability of an object at a given redshift being lensed, which we can then compare to optical 
depths for quasar-galaxy lensing, which has already been observed for about a dozen systems and should occur on 
the order of hundreds of cases in the SDSS. 

To estimate the string lensing probability, we require two pieces of information: the projected angular density of 
string and a lensing cross section for that string. With regards to the former, numerical simulations of string networks 
in an expanding universe can give reliable estimates of the string density pi s in a horizon volume, characterized by 
the horizon radius dn- Using the results of Bennett and Bouchet Q as a representative example, we see that in the 
matter dominated epoch the length of string in horizon units is a constant given by 

Lu = P_uA_ = z l± ^ (1) 

where fi is the energy per unit length of string. We shall treat this string as an ensemble of small links, with length 
Li in horizon units, distributed with a uniform probability density related to the above result, each with a random 
orientation with respect to the line of sight. These assumptions are appropriate when considering an average over 
many string networks. We shall also assume that each of the links is static, since we lack good information on the 
distribution of link velocities. In essence we are disregarding the effects of the Lorentz contraction (see next section 
for details), which means our estimate is probably only accurate to a factor of a few. We now subdivide space into 
differential volume elements such that the probability that more than one link, located by its center of mass, resides 
in the same volume is vanishingly small in comparision to the probability that one link resides in that volume. For 
the zth volume element, we find that the differential angle of sky subtended by the string is given by 

de i = n i ^-sin(0 i )dV i , (2) 

where u, is the number of links in the volume, f3i is a random orientation angle associated with each link, and d A is 
the angular diameter distance to the link 
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assuming a flat, matter dominated universe. Under the same assumptions, the volume clement and horizon distance 
can also be expressed in terms of the redshift: 
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To estimate the lensing cross section, we assume that each link has the same cross section as if it were part of an 
infinite straight string. In this case, the angular cross section per radian of string is 
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The distances D\ q and -D 9 are respectively the angular diameter distances from the lens to the quasar and from the 
source to the quasar. In flat space for a source located at a redshift of z 2 and an observer located at z%, the general 
form of the angular diameter distance is given by 
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We ignore the effects of structure formation on theses distances which introduce deviations from homogeneity that 
can effect the path length. Convolving the cross section with the angular string density, we can determine the optical 
depth for a quasar at a redshift of z — z q : 



(z q ) = 87rGfj,Li ^— Y] rii sin 2 (/?» ) d f D J. q 
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where Qo is the observed fraction of sky. Taking the expectation value for t, we find that (rii) = L; s / (Lidjj) and 
(sin 2 (/3)) =2/3. This leaves us with the integral over the observed volume 
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which can be performed analytically, giving the final result: 
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Now we would like to estimate the variance in r which will permit us to set limits on string parameters using lensing 
statistics. Results from numerical simulations indicate that on scales of about 0.1d#, the string undergoes a random 
walk. To make the problem tractible, we shall continue to treat each of these segments as uncorrelated with the rest, 
and presume that the random walk correlation is fairly small. The mean of the square of the optical depth is given by 
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The expectation value (riirij) in the limit of differential volumes is equal to 
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when rii and rij are uncorrelated, and n = (rii). Given this relation, one can easily show that 
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which has the analytic result 
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To get the full error in t we convolve this result with the theoretical uncertianty in the long string density ai s given 
in eq. (|l|), which adds in quadrature. Now consider a distribution of sources with mean number density as a function 
of redshift N(z). The expected number of observed lenses is given simply by 
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To calculate the variance, we must account for both the variance in r and the Poisson fluctuations in the source 
distribution. Including both these effects, we find the full variance is 
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a 2 = JdV (n(z)t(z) + N(z) 2 a 2 (z) + N(z) 2 T 2 (z)^j . (15) 

As a toy example, we consider a quasar distribution given by N(z) = S(z — 2)10 5 /7r, which rough approximates that 
of the Sloan Digital Sky Survey. The SDSS will observe one quarter of the full sky, so using the results of this section, 
we find for G[i — 10 -6 , that the number of observed quasar lenses is 18 ± 6. 

In figure [l] we show the optical depth for long cosmic strings, assuming a value Gfi = 10~ 6 , compared with the 
optical depth resulting from galaxy lenses as calculated by Turner, Ostriker, and Gott [|J. The galaxy estimate is 
based on an isothermal sphere model, which represent an upper limit to the lensing probability, as including a finite 
core to the galaxy tends to reduce the optical depth by up to 50% ||. Thus we expect that for a wide field survey 
anywhere from 10% to 30% of the observed lenses will be the result of long cosmic strings when compared to galaxy 
lenses (assuming G/i = 10~ 6 ). Given that order 10 quasar lensed have been observed, one may expect that a few 
string lenses should have been seen. In fact there are no lenses that have yet been ascribed to cosmic strings, but this 
is statistically unsuprising because of the large variance associated with small number statistics. We can, however, 
reasonably conclude that string tensions G/i greater that a few 10 -6 , which would predict tens of observed lenses, 
probably are ruled out, consistant with the limits from pulsar timing. 



III. GEODESIC DEFLECTION BY A LONG COSMIC STRING 

In a previous paper we derived an equation for the deflection of a null geodesic-corresponding to a photon 
trajectory in the geometric limit -arising from the gravitational field of a cosmic string loop. For long cosmic strings, 
we can make use of some of that derivation, but we must now be careful, as the long strings stretch to horizon scales, 
requiring us to consider certain surface terms which could safely be ignored when examining small loops. Let us begin 
again by assuming a weak field so that the full metric may be expressed as 

g l iu = n fJ ,u + h ll , v , (16) 

where r]^ — diag(— 1, 1,1,1) is the usual Minkowski metric and /i M „ is a small perturbation such that all terms of 
0(h 2 ) are negligible. For simplicity, we choose to work in the harmonic gauge, which implies the condition g^T^ = 0. 
Using this gauge choice, to linear order we are left with a simple wave equation for the metric 

□ 2 V = -IQttGS^, (17) 

where — — l/2rj^ u T^ is the traceless component of the stress energy tensor T^. If we decompose the photon 
four velocity 7 M into its zeroth and first order pieces, 7^ and 7^ respectively, then it is a straightforward calculation 
to solve the geodesic equation and show that the first order deflection for a photon emitted at ti and observed at t 2 
is given by 



12 



7a(l) = o / at ^"7(0)7(0) ~ V*7(o) 
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(18) 



Where we are integrating over the photons zeroth order trajectory, i.e. x^ = x^q + 7^(0) £■ For loops, which are 
compact, we could safely ignore the second surface term, but for now we should retain it when considering long 
cosmic strings as they are not compact objects. 

To make further progress, we need to consider the form of the stress energy tensor resulting from a cosmic string. 
Strings are well approximated as lineal gravitational sources, so the string configuration at any time is given by a two 
parameter vector function f(a,t) where t is the time and a is a parameter which runs along the conformal length of 
the string. One may straight forwardly infer from this that the string traceless stress energy may be written in the 
form 



J da F^\x- f(a,t)), (19) 



Where fi is again the string energy density and F^ v can be expressed in terms of / and its derivatives. Now, let us 
designate 7 a to be the contribution to the first order deflection which comes from the integral in eq. ([l§j), separating 
it from the surface term. Contracting this with a derivative, we get 
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= -8ttG/ S M „ 7 (Vro)> (20) 
where the second line comes from the metric equation ([l7|). Plugging in the stress energy from eq. (19), we see that 

d a ^ a = -SirGv [da ['dt F^5W(x +~f {0) t - f(a,t))rf 0)1 { 0) . (21) 



We can evaluate the time integral if we decompose / into components which are perpendicular, and parallel, /y, 
to the zeroth order photon trajectory, with the result 



d a j a = -8nGn / da 
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(22) 
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where to is the solution to the equation to = /||(c, to) — x\\q, i.e. the light cone time slice, and dots refer to derivatives 
with respect to time. Note that the limits on the a integral are constrained to the regions where a solution with 
t\ < <o < i-2 exists. Now we shall write out the left hand side of the equality in terms of a parallel and perpendicular 
gradients, d a j a ^ = Vi • / 7_L(i) +^|| '7||(i)> wnere the parallel gradient can be written as V|| •T'ii(i) = 75(o)7(o)<9' 5 7/3(i)- 
Now we consider the contraction ; 



7(0)7/3(1), wh ich from eq. (|18|) is 
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We are able to perform this integration because ifpftp is equivalant to taking a complete derivative with time, dj dt. 
Using the above result in conjunction with eq. (E2), it is easy to show 



Vj_ • 7j_(i) = — 87rG/i / da 



^7(0)7(0) 



5 {2 \x ±Q -fA 



2~dt 



/W(o)7(°) 



J t=t 



1 2 



(24) 



We can solve this equation by assuming that the first order deflection can be written as a gradient of a potential, 
7_l(i) = ^x^i leaving a two dimensional Poisson equation from which $ may be found by integrating over the Greens 
function for the two dimensional Laplacian, G(x±,x' ± ) = — hi(\x' ±Q — xj_ci\ 2 ) /An. Specifically, we get 
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Finally, to recover the perpendicular deflection, we take the gradient and add the surface term from eq. (|T^) which 
leaves us with 
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where h± is defined to be the perpendicular part of fyuaTrnv This result gives us exactly what we want for lensing 
calculations, the photon deflection away from its zeroth order path. The values of 70 and 711 which are equivalent to 
first order, give us the redshift of a photon as it passes a string, but we are not interested in this calculation here. 

Eq. (26) may at first seem to be a retrograde step as it involves a two dimensional integral of the metric in space 
where originally we had only a one dimensional integral of the metric over time. However, we argue that the first 
and third terms which contain the metric explicitly can be neglected. To do so we require an explicit solution for the 
metric in terms of the stress energy: 



(27) 

where t = t — \x — x'\ is the retarded time, and this solution is derived from the Greens function for the D 2 operator. 
Using our string stress energy given in eq. (19), we can reduce the metric to 
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Now considerpoints which are far from the string. The metric goes like an integral over |a; — /l^ 1 while the middle 
term in eq. (|26|) goes like an integral over \x±o — The latter represents the minimum distance between the 

zeroth order photon trajectory and the cosmic string, while the former will, in general, go like the distance from the 
photon to the string at the current time. In the case of string loops, we could expect that the two metric terms would 
fall off like the inverse of the distance while the middle term remained constant, thus allowing us to drop the explicit 
metric terms. With an infinite string, one must be more careful because the distance from the string can never be 
large with respect to the string size. However, when determining the image structure in lensing, it is not the absolute 
deflection of the photons which matters, but rather, it is the difference in deflection between two nearby rays that 
counts. In this case, we see that for photons far from the string, the difference in the contribution between the two 
metric terms declines with distance while the difference between the static terms remains constant. So, for photons 
which pass nearby the string, that is a small distance when compared to the source and observer, the effect of the 
metric terms then is merely to cause all of the images to be displaced, but not to alter the shape or relative orientation 
of the images. Thus, for the purposes of determining the structure of strong lensing due to strings, we may drop the 
explicit metric terms all together and write 
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where the omitted terms only add a constant to this result when the photon impact parameter is small compared to 
the distances of the source and observer to the lens. 



IV. BASIC GRAVITATIONAL LENSING THEORY 



Before we consider the specific example of a long cosmic string lens, let us first take a detour into the basic theory 
of gravitational lensing. For most calculations, it is more than adequate to consider photon trajectories as the path 
taken by null geodesies like those discussed in the previous section, so that the lensing system is described by simple 
geometry and wave properties are ignored. In figure ^ we show a pictorial representation of a gravitational lensing 
system with an optical axis defined to roughly intersect the center of the lens. A photon is emitted from a source 
S displaced from the optical axis by a vector rj at a distance of D s from an observer and Di s from the lens, and it 
travels in a straight line until it reaches the plane of the lens. There, its trajectory is deflected by a vector a, defined 
as the initial trajectory minus the final, and it again travels in a straight line to the observer located at a distance 
Di from the lens. This approximation is known as the thin lens approximation because the lens and the deflection it 
induces are presumed to occur in a single plane. To first order, this is valid for long cosmic strings. The location of 
the image /, where the ray intersects the lens plane, defines the vector £ which is measured from the the optical axis 
to / in the lens plane. Since we will be considering sources and observers that are separated on cosmological scales, 
The distances D;, Di s and D s are angular diameter distances (see eq. (|J)). 

In the small angle limit, the condition that a ray emitted from the source will have an image seen at / by an 
observer at O requires that 

v = ^a-D l Md). (30) 

J->1 
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This equation, often referred to as the lens equation, gives one the unique source location for any observed image, but 
its inverse, however, is not unique, as a single source may generate more that one image. It is convenient to recast 
the lens equation into a dimensionless form by dividing through by the length Di, so that the coordinates are given 
units of angular displacement with respect to the optical axis. Thus we define a new set of variables 

x = jr, (31) 

V_ 

D a ' 

_ _ Dis 
<* = <*— , 

(32) 

which yield a dimensionless lens equation 

y = x — ot(x). (33) 

The above equation can be inverted to give the location of images for a particular source location, but one can 
also use the information contained therein to determine the magnification of those images. Suppose a source emits a 
narrow pencil beam of photons which subtends a solid angle dQ* , while the image of the of the source beam subtends 
an angle dfi, so that the ratio of the solid angles dfl* /dfl will give the flux magnification of the image to the source. 
Using the lens equation, one can derive the magnification by considering the Jacobian matrix 

*> = & <"> 

and observing that the magnification factor Ai(x) = dQ*(x)/dTl(x) is given by the inverse of the determinant of Aij, 

M ^ = dctW < 35 > 

Note that this same factor will give the angular magnification for extended objects, so gravitational lenses conserve 
surface brightness. 

V. NUMERICAL CALCULATIONS WITH LONG STRINGS 
A. Constructing Long String Networks 



In § [II , we derived an expression for the photon deflection that will be relevant to gravitational lensing with cosmic 
strings; now we must consider the problem of generating realistic long cosmic strings. Analytic work, in particular 
scaling soultions, has successufully described some of the average properties of long strings, but it cannot give the 
details of the structure of a particular string. Only by simulating networks of interacting strings can we hope to 
generate realistic structure. This exercise can be greatly simplified if we restrict ourselves to flat, Minkowski space 
rather than considering the expanding universe, but the resulting structure of these strings will not quantitatively 
agree with those produced by expanding universe models. However, there should be good qualitative agreement, 
sufficient for our lensing analysis. That is, we expect that the structure of the lensed images should be similar enough 
to those which would result from more accurate string simulations to make general conclusions about string lenses. 

Let us begin by considering the equations of motion for a cosmic string. The dynamics of strings are dominated 
by their tension, while gravitational effects are suppressed by factors of order G(i, so for now we can ignore the back 
reaction, but we will discuss its effects later. We need the string location which is described by a four vector with 
a time like component f° =t and spatial displacement f(cr, t) that we used in the previous section. Then, considering 
only tension, the equations of motion for the string are 

F - = 0, (36) 

where dots refer to derivatives with respect to time and primes refer to derivatives with respect to a. Our choice of 
the harmonic gauge enforces two constraints, 
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f% = (37) 

and 

f 2 + f' 2 = 0, (38) 

which restrict the motion to transverse directions and mandate conservation of energy respectively. The stress energy 
can also be expressed as a function of / and has the form 

V =M / da (fj v - fX)5 (3) (x - /M))> (39) 

consistent with the form suggested in eq. (|l9|). 

For non-interacting strings, it would be sufhcient to specify some initial conditions consistent with the gauge 
constraints and evolve them using the wave equation. However, real strings can interact when two different segments 
intersect, causing them to reconnect with the opposite segment. To handle both the evolution and intersection of the 
string network, we shall turn to a clever algorithm first proposed by Smith and Villenkin jT). The foundation of the 
Smith Villenkin algorithm is the fact that for a set of points equally spaced in <r, separated by 5, the wave equation 
( |36j ) can be reduced exactly to a finite difference equation on a lattice of a and t points. In terms of the displacement 
vector /, we get 

f(a, t + S) = f(a + 6, t) + f(a - 5,t) - f(a, t - 6) (40) 
This second order equation can be reduced to a pair of first order equations if we consider the velocity, defined as 

f={f(a,t + 5)-l[f(a + 5,t) + f(a-6,t)]\/6. (41) 
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Thus, we get 



and 



f{cx, t + S) = ^ [f(<r + 6,t) + f(a - 5, t)] + f(a, (42) 

f(a, t + S) = \ \f(? + 6, t) + f(a - 5, *)] + [f(a + 26, t) - 2f(a, t) + f(a - 25, t)]/4S. (43) 

A complete solution can be specified by initially fixing the positions for every even point on the a lattice and velocities 
for every odd point. After the first time step, eq.'s ([42] & [43|) will give the positions for each odd point on the a lattice 
and velocities for each even point on the lattice. After the second time step, even points will again be positions and 
odd points will again be velocities, and so on, so that in general, the plane of a and t will be filled with interlocking 
diamond lattices of positions and velocities. 

The next challenge is to satisfy the gauge constraints on the lattice. In eq. (^) we have given a discrete velocity, 
and now we need a discrete version for the a derivative. It has the obvious form 



f'(a, t) = If (a + S,t)- f(a - S, t)}/26, (44) 
where the gauge constraints for the discrete / remain unchanged, 

/ • /' = 0, (45) 

and 

f + f 2 = l. (46) 

To ensure that these conditions would be preserved through the entire evolution of the string, Smith and Villenkin 
proposed discritizing space in the same way as a and t, that is, the space itself is a lattice of points with spacing 
8. String configurations are then described by a series of connected links between successive string points on the a 
lattice consisting of three possible types. The first type is a static link with f = and |A/ = 2|, implying that one 
coordinate of A/ is ±26 while the others are zero. The second type is a moving link for which |/'| = |/| = y/2/2. 
Two components of A/ are ±6 while the final is zero, and the velocity / also has two non zero components which 
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are each ±1/2 and must be normal to A/. The last type of link is a cusp where A/ = 0. The velocity of a cusp is 
1 corresponding to a link traveling at the speed of light parallel to one of the axes. One can verify easily that all of 
these links satisfy the gauge conditions and that the constraints will be preserved by the equations of motion. 

Finally, to accurately evolve a string network, we must account for string inter-commutations which occur when 
different segments of the string collide. We have defined links which are connections of two successive string positions 
/(cr, t) and f(a + S, t). A collision occurs when two string points /(crj, t) and f(<jj,t) both fall on the same location. If 
this happens, we inter-commute the links so now the point / (<7»,t) is connected to f(<Jj +S,t) and the point f(aj,t) 
is connected to /(<7j + 5, t) with care taken to ensure that the proper velocity is assigned to each link. Prior to each 
time step, one tests each of the points to see if it lies in the same position as any other. One inter-commutes all 
the colliding links, and then evolves time for one step S and checks the strings again for inter-commutations. To 
ensure that cusps are not mistaken for collisions, we set a minimum separation in a space, 4(5, required before an 
inter-commutation on a distinct string is allowed. This fixes a minimum loop size, but the structure of the long strings 
is insensitive to the choice so long as it is small in comparison to the overall string length. In passing we mention 
that the order N 2 process of testing each point for intersection can be reduced to an order few N log N process if one 
sorts the points with respect to their positions using a quick sorting routine first. 

The Smith- Villenkin algorithm is a powerful tool for evolving networks of cosmic strings, and now we shall consider 
the initial conditions which we will use to produce long string segments for gravitational lensing. We start with the 
standard initial conditions introduced by Vachaspati and Villenkin who lay down a periodic lattice of phases, and 
locate strings through the faces of the lattice cubes which have non zero winding number. The result is a network 
of strings made of static segments which are parallel to one of the axes of the box, ideally suited as initial conditions 
for the Smith- Villenkin evolution algorithm. To ensure that space discritization effects do not strongly influence the 
resulting strings, each of these segments was subdivided into sixteen static links 2d long, where the overall segment 
length is 32(5. Tests of our simulations found that subdividing the segments with a larger number of links did not 
change structure of the evolved network, so spatial discritization effects have been minimized. Ideally, we would like 
to evolve this network until it has completely relaxed, and use the resulting long strings for our lensing calculations. 
Unfortunately, the periodic boundary conditions insure that there can be no net string flux through the box, meaning 
that given sufficient time, all of the long strings will fragment into loops leaving us nothing to study. To minimize 
these periodic effects, we evolve the network for a time equal to half the light crossing, i.e. if the box is n5 wide, we 
take n/2 time steps, each <5 long. We expect scales on the box size to preserve the structure of the initial conditions, 
but on scales smaller than nS/2, interactions should have sufficient time to relax the system to its final state. 

Previous analysis of evolved string networks has shown that the small scale structure of the surviving long strings 
is self-similar, that is, the structure is statistically the same on all scales well below the horizon We compare 
our results with our predecessors, as a consistency check. The simplest test is to measure the fractal dimension of 
the string defined as the exponent n such that L cx D n , where L is the mean conformal length of string measured 
between points separated by a distance D in physical space. For self-similar structures, n is a constant for all D, and 
for the particular example of a random walk string, n = 2.0. In figure || we show a log-log graph of the conformal 
length L plotted as a function of D for the long strings produced in a box simulation with periodic length 1024(5 at 
times t = OS, 256(5, and 512(5. Note that our initial segments were formed in a 32 3 box of phases and each initial 
segment was 32(5. We can see that for the unevolved strings that the fractal dimension is roughly uniform above the 
initial link scale. A linear regression of these points gives a slope of n = 2.0 which is what one expects from a random 
walk. As the simulation evolves, we see that the shorter length scales begin to relax into a new structure with a 
smaller fractal dimension, and by the time t — 5125, scale below D ~ 512(5 have almost completely relaxed. A fit to 
these points gives a fractal dimension of n = 1.3, a result consistent with those of Sakellariadou and Villenkin |jTo| , 
and interestingly, also consistent with the early time results of expanding universe simulations. In this last type of 
simulation, apparently the strings first fragment and relax on time scales that are short compared to the expansion 
time of the universes and then are stretched so that n falls below the initial flat space value. 

B. Finding Images With Long Strings 

Now that we have a network of realistic long strings to work with, we consider how to use them in gravitational 
lensing systems. The best approach would be to evolve a large network with very fine spatial resolution (on the 
order of 10~ 5 the box length, where we could associate that scale with a few times the horizon scale), and just send 
photons through the box. Of course, the numerical resolution required to perform such a calculation is well beyond 
the capacity of current computers, so as an alternative, we exploit the fractal nature of the small scale structure of 
the string discussed in the previous subsection. Sections of string which can fit in a box of length about half the 
simulation box relax into a self similar structure described by a constant fractal dimension, meaning that if we were 
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to magnify a small piece of this string, we would observe the same structure in proportion to the new scale. Thus, 
to use our long strings which are not resolved on the scales relevant to gravitational lensing, we need only rescale the 
string, as long as we restrict ourselves to segments which have had sufficient time to relax. 

Let us be more specific about precisely how we accomplish this. From our simulations described in the previous 
section, we have a periodic box filled with loops and long strings. We can eliminate the loops by considering only 
strings which are significantly longer than the box length, leaving the long strings which tend to wrap around the 
box periodically once or more. It is the smaller scales, for which the string has relaxed to its final structure, that 
we wish to consider. To do so, we remove a long string from the box entirely, laying it end to end no longer in 
a periodic box. If the string wraps around the original network box more than once, we connect the the box size 
segments by shifting the endpoints to make one super-long string. Now we need to determine the the light cone 
projection of the string required to find the deflection given by eq. (|2^). We accomplish this by allowing our single 
long string to evolve independently of the network, turning off all inter-commutations and connecting the end points 
periodically. Inter-commutations are ignored because we do not want to alter the small structure of the string. A 
photon is presumed to travel along one of the axes, and the location of the intersection of each of the string points 
with the light cone, along with the string velocity is recorded. Since each string point is equal spaced in a, we can use 
the set of light cone projected points to reduce the integral in eq. ( p6| ) to a discrete sum. We do not, however, wish to 
use the entire long string, since the structure has not relaxed on the largest scales. Instead, starting at an arbitrary 
point, we select shorter segments-those which can fit in a box with side length half that of the original network-and 
use only these points to calculate the photon deflection. 

In truncating the summation and considering only a finite string segment, we obviously introduce errors which we 
would like to quantify. As a order of magnitude estimate, let us consider the special case of an finite straight string 
segment perpendicular to the photon trajectory with equal lengths I above and below the photon axis. We consider 
the photon to be moving along the z axis and place the string segment parallel to the y axis a distance d from the 
origin (defined by the photon) so that the deflection will be along the x direction. Using eq. (^9|), while observing 
that /|| = and F^j^j" = 1, we see that the magnitude of the deflection will be proportional to 

J\ ^T^2 =2tan-V/d), (47) 

and the truncation error, that is the difference between I — > oo and finite £, goes as it — 2ta,n~ 1 (£/d). In the limit of 
large I this is approximately 2d/£. If instead we wish to look at a fractal string, then we should replace the a 2 in the 
integral with something proportional to cr 2 /™ which gives a truncation error which falls like (g?/£)( 2 /") -1 . So, we see 
that the error depends on how far the photon passes from the string, which is fortuitous since most of the interesting 
lensing occurs for photons which pass close to the string segment. We can also see that the thin lens approximation 
is justified. So long as d is much smaller than the distances of the source and observer from the lens,~ D, then the 
deflection will occur in a region within a few lengths d from the string which, and including the effects arising from 
the non-instantaneous deflection would be second order in d/D. 

Given the projected string segment, we have everything necessary to examine gravitational lensing, but now the 
challenge is to solve the lens equation. Currently there is no fast way to invert a multidimensional equation like eq. 
(|33|), so our only choice is to solve the problem by brute force. Given a source and lens redshift, we calculate the 
deflection a(x) on a uniform grid in x space, and then use the lens equation to solve for y at each x. In other words, 
we are mapping a uniform grid in the image plane back onto the source plane where we can use this map to locate 
the images of a particular source. We compare triangles made of nearest neighbors in the image grid mapped onto 
the source plane to the locations of the sources. If a source point is enclosed by a mapped image triangle, then its 
image must lie somewhere in that triangle in the image plane, and thus, for any source point, we can locate its image 
to the uncertainty given by the image grid spacing. One must be careful, though, when considering image triangles 
that enclose a piece of string because the photon deflection is discontinuous as x crosses the string. Images in these 
triangles should be ignored because they represent photons which must pass through the string itself to be observed 
and are therefore spurious. Since our strings are resolved down to the scale 5, it is natural that the image grid spacing 
should be 8 as well. And, having assumed a self-similar structure for the string, we are free to choose the physical 
scale that 6 represents. For the objects we shall consider for lensing, we find that a good choice for this scale is 
S/Di — 0.1 arc sec. 

Using the techniques described above, we shall consider two types of sources, point -like and extended, which will 
roughly correspond to quasars and galaxies respectively. Let us first consider quasars as they represent what are likely 
the best objects to look at when trying to observe strings through lensing. Quasars possess a trio of virtues regarding 
lensing, namely they are high redshift objects (z ~ 1 — 5), they are bright and therefore easily observed, and they are 
typically separated by large enough distances that the chance coincidence of two different quasars being separated at 
the same scales as typical lensing systems is rare. However, there are a sufficient number to make the observation 
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of string quasar lensing systems probable (see iQ). Quasars are compact, distant, and cannot be spatially resolved 
with current technology. Thus we shall treat them as idealized point sources. The locations of the resulting images 
are found as we have described, along with their magnifications which are determined by eq. ([35j). Figures |[ ||, and 
H show the results quasar lensing with three different segments of cosmic string. In each panel, the source location is 
shown as a hatched circle, the images as open circles, and the projected string segment as a dashed line. The relative 
area of the image circles gives the ratios of magnifications of the image to the source. The string is located at a 
redshift of one and the quasar is located at a redshift of two, while the angular size of each panel is 25 arc sec. The 
full length of string used in the lensing calculations is not shown, but we use the same segments for the galaxy lenses, 
and we do show the full segment in figures ^ ^, and [ll] respectively. Typical in many of these examples are a number 
of small demagnified images which reside close to the string itself. These results are qualitatively similar to those one 
would see if we replaced the string by several point masses with masses similar to the energy in the string. With 
point masses, The deflection induced by neighboring points cancels around the midpoint between the two masses. 
Small images tend to form here because the deflection angle changes rapidly around the minima. For strings, it is 
the kinks and wiggles which provide a similar opportunity. Inside a kink, contributions from different parts of the 
string can cancel, leading to rapid changes in the deflection. Kinks can also produce small images by just being large 
concentrations of energy. They produce images in a manner similar to point masses when the source is outside the 
Einstein ring. The secondary image is smaller and therefore demagnified. 

High redshift galaxies also provide interesting candidates for cosmic string lensing, but observing these systems is 
significantly more challenging. The greatest difficulty lies in observing such objects because they are so faint. Things 
arc further complicated because foreground galaxies may also contaminate the systems making it more difficult to 
observe the images. Also, typical galaxies are not spherical making it difficult to determine if one is observing a 
lensing arc or merely an edge on galaxy. However, to qualitatively illustrate a galaxy lensing system, we consider 
an idealized case of a set of extended spherical sources located at a redshift of z = 2, lensed by a string at z = 1. 
The angular size of the sources is 0.5 arc sec which corresponds roughly to the observed size of high redshift galaxies. 
Seven such sources are scattered randomly in side a (25 arc sec) 2 viewing area mimicking the approximate angular 
density of real high redshift galaxies j[l) . In figures 0, and [n] we show the full string segment used in the lensing 
calculation along with the sources to be lensed. The dashed box shows the area for which we calculate the images, 
where outside this region, we expect that the photon deflections will not be particularly accurate due to truncation 
error. In figures g [l0| and 12 we show the observed images corresponding their respective sources, and again we note 
the proliferation of smaller, demagnified images similar to those observed in the quasar images. 



VI. CONCLUSIONS 



From the results in the previous section, we can see that cosmic strings can produce images which have characteristics 
unlike the more prosaic galaxy lens. In particular, the proliferation of small demagnified images is a signature that 
may be unique to long cosmic strings. This suggest the exciting possibility that a cosmic string can be positively 
identified through gravitational lensing, confirming the existence of these topological defects. Unfortunately, our 
enthusiasm must be tempered by two important caveats. The first is that these strings are the product of Minkowski 
space simulations and do not include the effects of the expanding universe. When theses are considered, the structure 
on the string is stretched so that the fractal dimension falls to approximately n = 1.1 Jl2), so we expect real strings 
to have smaller kinks and wiggles. The other effect that we ignored is that of gravitational radiation. For loops it is 
well known that they radiate energy at a rate of TG/i where T is on the order of 100, so loops shorter than TG/xt will 
have radiated completely away. The effects of gravitational rediation on the structure of fractal strings is not as well 
understood. An alalytic example of a helical string has been calculated by Sakellariadou |ll| while the case of small 
amplitude kinks has been considered by Hindmarsh |l4|] . Their results suggest that fluctuations on long strings should 
have a life expectancy of about d/Gfi where d is the typical separation length between kinks-on the same order as the 
fluctuation. So, one might expect that the strings are straight on scales smaller than Gfj, times the horizon scale at 
the epoch of lensing. For a string located at z = 1 this corresponds to an angular scale of 0.5 arc sec, so gravitational 
radiation may just influence the long string lensing structure. Our examples then represent the most extreme results 
that one should expect, and real signatures may be less distinct. However, because the small scale structure of our 
strings is responsible for the demagnified images, we can conclude that gravitational lensing may be a good way to 
probe that small scale structure. If strings are relatively smooth on the scales we considered in our examples, then 
any strong lenses observed should produce a pair of undistorted images like those resulting from an infinite straight 
string. Conversely, if there is significant small scale structure, then one expects to see a number of demagnified images 
like those seen in the figures. 

In §JllL we found that quasar-string lensing systems should constitute about 10% to 30% of the observed quasar 
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lenses, where the rest arise from lensing with galaxies. We suggest that the search for gravitational lenses could on 
its own confirm or rule out the cosmic string model of structure formation depending on whether such lenses are 
observed. Since the strings are obviously correlated, a large angle of sky coverage with a large number of quasars is 
required. In fact, precisely such a survey is currently being developed, namely the Sloan Digital Sky Survey. Some 10 5 
quasars with a mean redshift of two are to be observed in a tt steradian slice of the sky. We have estimated that one 
should observe on the order of ten string lenses in the SDSS for an f2 = 1, cosmic string model. The failure to observe 
any string lensed quasars would require G[i < 10 6 , making strings an unlikely candidate for structure formation. 
Conversely, should a suspected string lens be observed, it is possible to confirm it by looking at galaxy observations 
in the neighborhood of the lens. With precise observations, one expects to see lensing of high redshift galaxies from 
nearby parts of the same string. In concert then, quasar observations followed by galaxy observations could provide 
definitive proof of GUT scale cosmic strings, or rule out cosmic strings as a viable model for the formation of large 
scale structure. 

We would like to thank Tanmay Vachaspati for his help and suggestions. This work was supported with funding 
from the Department of Energy. 
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FIG. 1. The probability of a source being lensed by a long cosmic string plotted as a function of source redshift (solid). For 
comparison, the probability of a source being lensed by a galaxy is shown (dashed). 
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FIG. 2. A schematic representation of a gravitational lensing system. 
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FIG. 3. String conformal length L plotted as a function of point separation distance D. The solid curve is for t = 0, the 
dotted for t = 2565 and the dashed for t — 5125. The periodic length of the box is 10245. 
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FIG. 4. A string-quasar lensing system. The string is shown as a dashed line while the source is shown as a hatched circle 
with images shown as open circles. The relative areas are indicative of the magnifications of the images. 
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FIG. 5. Same as Fig. 
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FIG. 6. Same as Fig. 
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FIG. 7. A galaxy-string lensing system. Shown is the full string used in calculating the photon deflection along with several 
galactic sources. The dashed box shows the region in which the images of these sources are calculated. The string is located at 
z = 1 and the galaxies at 2 = 2. 
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FIG. 8. The images observed for the lensing system shown in Fig. [?] 
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FIG. 9. Another string-galaxy system like Fig. [?] 
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FIG. 10. The images observed for the lensing system shown in Fig. ^| 
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FIG. 11. Another string-galaxy system like Fig. ^ 
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FIG. 12. The images observed for the lensing system shown in Fig. [ll] 
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